Evolution of the Dark Matter Distribution at the Galactic Center 
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Annihilation radiation from neutralino dark matter at the Galactic center (GC) would be greatly 
enhanced if the dark matter were strongly clustered around the supermassive black hole (SBH). 
The existence of a dark-matter "spike" is made plausible by the observed, steeply-rising stellar 
density near the GC SBH. Here the time-dependent equations describing gravitational interaction 
of the dark matter with the stars are solved. Scattering of dark matter particles by stars would 
substantially lower the dark matter density near the GC SBH over fO Gyr, due both to kinetic 
heating, and to capture of dark matter particles by the SBH. This evolution implies a decrease by 
several orders of magnitude in the observable flux of annihilation products compared with models 
that associate a steep dark matter spike with the SBH. 



Neutralinos in supersymmetry are likely candidates for 
the non-baryonic dark matter |lj, |2( • If neutralinos make 
up a large fraction of the dark matter in the galactic 
halo, pair annihilations will produce an excess of pho- 
tons which may be observed in gamma ray detectors Q . 
The galactic center (GC) is a promising target for such 
searches since the dark matter density is predicted to 
rise as p ~ r -1 at the centers of dark matter halos Q. 
In addition, the GC contains a supermassive black hole 
(SBH) with mass M. w 1O 6 5 M 0. "Adiabatic growth" 
models in which the SBH remains stationary as it grows 
predict the formation of a steep power-law density pro- 
file around the SBH, a "spike," and an increase by many 
orders of magnitude in the amplitude of the neutralino 
annihilation signal Q. 

The bulges of galaxies like the Milky Way are believed 
to have formed via mergers of pre-existing stellar sys- 
tems. If the latter contained SBHs, a merger would result 
in the formation of a binary SBH |jj . The density of stars 
and dark matter around a binary SBH drops rapidly as 
the binary ejects matter via the gravitational slingshot 
||. Even a binary with mass ratio as extreme as 1 : 10 
would efficiently destroy a dark matter spike on parsec 
scales 0. 

Evidence for the scouring effect of binarySBHs is seen 
at the centers of the brightest galaxies [lJJ,lll|, where the 
stellar density profiles are nearly flat and sometimes even 
exhibit a central minimum |12j . However in fainter ellip- 
tical galaxies and in the bulges of spiral galaxies like the 
Milky Way, steeply-rising stellar densities are observed: 
~ r~ 7 ,1.5 < 7 < 2.5. In these galaxies, the most 
recent mergers may have taken place before the era at 
which SBHs formed, allowing the stellar density near the 
SBH to remain high. Since stars and dark matter re- 
spond similarly to the presence of a SBH, galaxies with 
steeply-rising stellar densities are the most plausible sites 
for steeply-rising dark matter densities and hence for the 
detection of annihilation radiation [l3j | . 

Stars near the SBH would also interact with the dark 
matter via gravitational scattering 0, [^. Here, the 
time-dependent equations describing the scattering of 
dark matter particles off of stars in the presence of a SBH 



are solved. Scattering decreases the density of a dark 
matter spike by kinetic heating, and by driving particles 
into the SBH. The result, after 10 Gyr, is the virtual 
dissolution of the dark matter spike in a galaxy like the 
Milky Way. This result suggests that enhancements in 
the dark matter density around the GC SBH would be 
modest whether or not the Milky Way bulge has experi- 
enced the scouring effects of a binary SBH. 

Let rn be the radius of gravitational influence of the 
SBH, with M*(r < r h ) = 2M.. For the GC SBH, r h « 
1.67 pc [lfj. After growth of the SBH (assumed to remain 
fixed with respect to the bulge) , the dark matter density 
is approximately 



p(r) — p(r b ) x (r/r b ) 7sp , r < r b 
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where 7 sp = 2 + 1/(4 — j c ), poc r 1c - is the dark matter 
density before growth of the SBH, and r b « 0.2r^ [l7| . 

An estimate of the local heating rate of dark matter 
particles due to gravitational encounters with stars is the 
change per unit time of e = ^mv 2 ms , the mean kinetic en- 
ergy of the dark matter particles. Assuming Maxwellian 
velocity distributions for the stars and dark matter, 
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where e* = \m i ,v 2 krms , p± is the stellar mass density 
and In A is the Coulomb logarithm [l^. Taking the 
limit m <C m+ and assuming v rms ~ D^ ms , appropri- 
ate shortly after the dark matter spike forms, the local 
heating time becomes 
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the latter expression uses the observed stellar mass den- 
sity near the GC SBH, p+ » 3.2 x lO 5 Af pc- 3 (r/lpc)- 7 , 
7 = 1.4 ±0.1 16j , and v+, rms « 1.12(GM./r) 1 / 2 with 
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M. = 3 x 10 6 M Q . The Coulomb logarithm was set to 
In A = ]n(0AN) with N w 6 x 10 6 the number of stars 
within jl9j . The time to heat the dark matter is nearly 
independent of radius and shorter by a factor ~ 5 than 
the age (~ 10 Gyr; |2j}) of the stellar bulge. 

The change with time of the dark matter density can 
be computed from the Fokker-Planck equation describ- 
ing the evolution of f(r,v,t), the mass density of dark 
matter particles in phase space, due to gravitational in- 
teractions with stars [21| . We assume that / is isotropic 
in velocity space, / = f(E,t), with E = — v 2 /2 + 4>(r) the 
binding energy per unit mass of a dark matter particle, 
and <p(r) = —$(?•), with <E>(r) the gravitational potential 
due to the SBH and the stars. The kinetic equation de- 
scribing the evolution of f(E,t) due to scattering off of 
stars with masses 3> to is 
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Here p{E) = Ay/2 J^ {E) drr 2 ^/ <j>(r) — E = -dq/dE 
is the phase space volume accessible per unit of energy, 
with 4>(r max ) = E. The mass density of dark mat- 
ter particles in £-space is N(E)dE = An 2 p(E)f(E)dE. 
The dark matter heating rate is determined by h+^E) = 
J m*f+(E, to*)c£to* with f+(E, m*)dTO* the mass density 
of stars in phase space in the interval to* to to* + g?to*. If 
the distribution n(TO*)d?7i* of stellar masses is assumed 
independent of energy (i.e. distance from the black hole), 
then h+(E) = to*/*(.E), with /*(-E) the total mass den- 
sity of stars in phase space and 
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Fi c (E) is the flux of stars that are scattered from low 
angular momentum orbits into the SBH and is discussed 
in more detail below. Eqs. (0} assume that small-angle 
scatterings dominate the evolution of / and that the grav- 
itational potential changes on a time scale long compared 
with Xiocai; the latter assumption is valid as long as the 
gravitational acceleration is not produced dominantly by 
the dark matter particles themselves. 

Neglecting F lc , the total energy £ = /jf 2 N(E)EdE of 
dark matter particles in the energy range E\ < E < E% 



changes with time according to 
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Q(E) = 167r 2 G 2 TO*lnA / dE'f*(E'). (6d) 
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The third term in eq. JHt); which is always negative, 
represents heating of the dark matter. We accordingly 
define the (non-local) time scale for heating of the dark 
matter particles to be 
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This expression may be used to estimate the dissolution 
time of a dark matter spike. Assume that both stars and 
dark matter particles initially have power-law density 
profiles near the SBH: p(r,t — 0) oc r~ l3V , p*(r) oc r~ 7 , 
r ^ r h, an d that <f)(r) = GM,/r. The isotropic dis- 
tribution function corresponding to an r ' density pro- 
file in an r" - 1 potential is f(E) oc E^-' i l 2 . Setting 
Ei = 4>{rh) = GM.Jrh and E 2 — > oo, T heat becomes 
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When 7 = 3/2, equal within the uncertainties with the 
slope of the stellar cusp around the Milky Way SBH ^(| , 
the coefficient A(^,j sp ) in eq. JSJ is independent of j sp 
and 
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The effective stellar mass to* that appears in equa- 
tion (9a) depends via equation © on the mass func- 
tion n(m*) of stars in the GC stellar cusp. While ti(to*) 
is not strongly constrained by observations, either at 
the hi gh or low mass ends, it is sometimes assumed 
(e.g. |l6j|) to be a power law with Salpeter in- 
dex, ti(to*) oc to~( 1+q ), a w 1.35. Setting TO* im i„ = 
0.08M© and Tn+ max = (5)1O(2O)M then yields to* w 
(0.8)1. 2(1.8)M . 
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FIG. 1: Evolution of the dark matter phase space density 
f(E) and mass density p(r) due to gravitational interactions 
with stars around the GC SBH. The initial dark matter den- 
sity is given by eq. 1 with 7 sp = 7/3 and 7 C = 1. Times 
shown are r = (heavy curves) and t = 2,4, 20 where r 
is the time in units of Theat (eq. [5J ; 10 < r < 20 corresponds 
roughly to the age of the galactic bulge. 



We take Th oa t as defined in eq. © as our unit of time in 
what follows, with t = £/Th ca t- The age of the majority 
of the stars near the GC is > 10 Gyr 20], although some 
much younger (t < 10 7 yr) stars are present in the cusp at 
distances < 0.1 pc from the SBH [11 Setting t = 10 
Gyr, r h = 1.67 pc and 1 < m*/M < 2 gives 10 < 
t < 20. If the young stellar population is continually 
replenished, m* and r could be larger, implying a higher 
mean rate of dark matter heating. 

Diffusion in energy will cause a modest loss of stars 
to the SBH, M = -F E (E 2 ), E 2 ~ c 2 . A much greater 
capture rate is implied by scattering of dark matter par- 
ticles on low angular momentum (eccentric) orbits into 
the SBH |25j. The loss rate is given approximately by 




F lc (E) « 4n 2 P(E)J 2 (E)Jl(E)R^ (10a) 
» S(E)f{E), (10b) 
S(E) = 4n 2 P{E)J 2 {E)jl{E)[\nR (E)- 1 ]^ 1 (10c) 

Here R = J 2 /J 2 (E) is a scaled angular momentum 
variable with J C {E) the angular momentum of a circu- 
lar orbit of energy E; P is the period of a radial or- 
bit; 71 is the orbit-averaged angular momentum diffusion 
coefficient ((AR) 2 ) /2R; and R is the value of the an- 
gular momentum variable at which / drops to zero due 
to the competing effects of capture and diffusion. The 
final line of eq. (| 1 Oil assumes / ~ \n[R/ Rq(E)] near 
the loss cone [2(| ]22|. Cohn & Kulsrud J23] give ex- 
pressions for i?o as a function of E based on solutions 
to the R— dependent Fokker-Planck equation; we adopt 
their expressions here. The angular momentum diffusion 
coefficients used by these authors, for modelling systems 
containing a single stellar mass, may be shown to remain 
unchanged when the scattered objects (here dark mat- 
ter) have masses much less that those of the scatterers 
(stars) . The rate of loss of stars predicted by eq. I|10|) de- 
pends only weakly on the radius of the capturing sphere, 
which we set to 2GM,/c 2 . The diffusion coefficient JI is 
of order T^ at , hence F lc (E) a N(E)/T hc , t {E) lni?^ 1 . 

The detailed evolution of the dark matter density 
around the GC SBH was computed by integrating eq. 
(0Ji) forward in time. The stellar density was modelled 
via Dehnen's [3] density law, p*(r) cx (?'/ro) _7 (l + 
r/ro)~ / ~ 4 with 7 = 1.4 and ro chosen to match the ob- 
served stellar density at r < rh 16] . Fig. 1 shows the 
evolution of the dark matter density assuming 7^ = 7/3 
and 7 C = 1, the values corresponding to a spike that 
developed in response to adiabatic growth of a SBH in 
a dark matter halo with j c = 1.0. Fig. 1 shows that 
scattering of dark matter particles by stars causes the 
dark matter density to drop and the spike to flatten, 
within a radius ~ rh where the heating time is shorter 
than the age of the bulge. Fig. 2 shows the dark mat- 
ter density at r = 10 for a range of initial spike profiles, 
1.0 < 7sp < 2.75. The mean density within rh drops 
by several orders of magnitude when 7 sp « 2. When 
ls P = 7c = 1, i-e. no spike, the net effect of the heating 
is to increase the dark matter density slightly. 

In the absence of scattering into the SBH, eqs. Pfl 
have the time- independent solution f(E) = constant, 
p ~ r _3//2 0] ■ As Figs. 1 and 2 show, there is a tendency 
to evolve toward this characteristic profile, although a 
number of factors keep it from being precisely reached, 
including the finite evolution time; the presence of the 
loss term Fi c ; and the fact that / = constant can only 
hold true over a finite range of energies given the bound- 
ary conditions on /. Nevertheless, at late times (r > 20), 
the solutions found here are generally well described by 
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FIG. 2: Dark matter density at r = (dotted lines) and 
t — 10 (solid lines) given an initial density that satisfies eq. 1, 
with 7c = 1 and 7sp = (1.00, 1.50, 1.75, 2^00, 2.25, 2.50, 2.75). 

p ~ r -3/2 at radii 10 -5 < r j Th < 10 -2_ 

The flux of dark matter annihilation photons along 
a direction that makes an angle tp with respect to 
the GC is proportional to the line-of-sight integral 
J^p 2 dl. Following earlier authors we define the di- 

mensionless form factor J(ip) = K J^p 2 (l)dl, K^ 1 — 

(8.5 kpc)(0.3 GeV/cm 3 ). Given a photon detector with 
angular acceptance Aft directed toward the GC, the sig- 
nal is proportional to 

(j) = ^J AQ m<M. (id 

Figure 3 shows the evolution of (J) for A£l = 10 -5 (1CT 3 ) 
sr; the first value is the approximate solid angle of the 
detectors in GLAST |29j and in atmospheric Cerenkov 
telescopes like VERITAS |3(J, while the larger angle 
corresponds approximately to EGRET |3J. The dark 
matter density was normalized to a fiducial value of 
p = 100A/ Q pc~ 3 at r = r/,; J and (J) scale as p 2 (r^). 
We note that p{rh) is uncertain and could be much lower 
0, H3, H^|. Figure 3 shows that the very large initial 
values of (J) are rapidly diminished as the spike is dis- 
solved; by r = 10, (J) has dropped below ~ (10 4 ,10 3 ), 
Aft = (10" 5 ,10- 3 ) for all 7sp < 2.5. These values are 
similar to what would be predicted for the central regions 
of a dark matter halo in the absence of a SBH [j| . 

The dark mass captured by the SBH after 10 Gyr is 
less than 10 4 M Q (p(r h ) = IOOMqPC" 3 ) for all the in- 
tegrations presented here. Various schemes have been 
discussed for increasing the captured mass in stars or 
dark matter to much greater values, perhaps of order 
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FIG. 3: Evolution of (J), a dimensionless measure of the cus- 
piness of the dark matter spike (eq. 1111 . for two values of 
the solid angle AQ of a detector centered on the SBH. Curves 
are shown for 7 sp = 1.50 (lower), 1.75,2.00,2.25,2.50, and 
2.75 (upper). The dark matter density was normalized to an 
initial value of IOOMqpc -3 at r = r^; (J) scales as p 2 (rh). 

M.. These include making the dark matter collisional 
|34| : assuming instantaneous replenishment of the loss 
cone [35l |3(| ; or allowing the stellar potential to be non- 
axisymmetric The first two mechanisms are ad hoc; 
the third, if it applies to the GC, might allow the persis- 
tence of a dark matter spike in the face of scattering and 
capture by increasing the mass in dark matter particles 
that can interact with the SBH. 
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